Gravitational collapse in asymptotically Anti-de Sitter/de Sitter backgrounds 
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We study here the gravitational collapse of a matter cloud with a non- vanishing tangential pressure 
in the presence of a non-zero cosmological term. Conditions for bounce and singularity formation 
are derived for the model. It is also shown that when the tangential pressures vanish, the bounce 
and singularity conditions reduce to that of the dust case studied earlier. The collapsing interior is 
matched with an exterior which is asymptotically de Sitter or anti de Sitter, depending on the sign 
of cosmological constant. The junction conditions for matching the cloud to exterior are specified. 
The effect of the cosmological term on apparent horizons is studied in some detail, and the nature of 
central singularity is analyzed. We also discuss here the visibility of the singularity and implications 
for the cosmic censorship conjecture. 
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I. INTRODUCTION 

The gravitational collapse of a matter cloud which is 
pressureless dust, and its dynamical evolution as gov- 
erned by the Einstein's equations, were first studied in 
detail by Oppenheimer and Snyder In recent years 
there have been extensive studies on gravitational col- 
lapse to examine the final fate of a collapsing cloud in 
order to investigate the final state of such a collapse in 
terms of formation of a black hole or naked singularity. 
These studies would throw important light on the nature 
of cosmic censorship and possible mathematical formu- 
lations for the same (for some recent reviews, see e.g. 
@). Such studies have already helped to rule out several 
possible versions of cosmic censorship, where a precise 
and well-defined formulation itself has been a major un- 
resolved problem so far. While understanding the nature 
of dynamical gravitational collapse within the framework 
of Einstein's gravity itself is a problem with considerable 
astrophysical significance, the understanding of cosmic 
censorship, if it is valid in some form, is another major 
motivation for such collapse studies. Most of these in- 
vestigations so far, however, have assumed a vanishing 
cosmological term (A) . 

The cosmological constant is sometimes thought of as 
a constant term in the Lagrangian density of general rela- 
tivity, and it is also theorized that A may be related to the 
energy density of vacuum (See Q and references therein). 
Recent astronomical observations of high red-shift type 
la Supernovae 0| strongly indicate that the universe may 
be undergoing an accelerated expansion and it is believed 
that this may be due to a non-vanishing positive cos- 
mological constant. On the other hand, the proposed 
AdS/CFT conjecture [5| in string theory has generated 
interest in the possibility of space-times with a negative 



'Electronic address: arunmadhav@gmail.com 
^Electronic address: goswami@tifr.res.in 
t Electronic address: psj@tifr.res.in 



cosmological constant. The conjecture relates string the- 
ory in a spacetime where the non-compact part is asymp- 
totically anti-de Sitter (AdS), to a conformal field theory 
in a space isomorphic to the boundary of AdS. Whereas 
the obvious effect of a positive cosmological constant is 
to slow down gravitational collapse, that of the negative 
cosmological constant is to supplement the gravitational 
forces. 

In this context, it seems pertinent to study the dynam- 
ical collapse of matter clouds when A ^ 0. Our purpose 
here is to study a class of collapsing models which incor- 
porate pressure, and which are asymptotically either de 
Sitter or anti-de Sitter geometry, depending on the sign 
of the cosmological term. Dust models with A are known 
in literature Q and there have been some studies on dust 
collapse with A in recent years (see e.g. Q)- We intend 
to discuss here a sufficiently general fluid model which 
allows pressure to be non-zero, and which also allows the 
cosmological term to be non- vanishing. Specifically, we 
study models with a non-vanishing tangential pressure 
pe |9, together with A ^ 0. Collapse models with a tan- 
gential pressure have been studied extensively, but not 
with a non-zero A §j. Allowing the collapse to develop 
from regular initial conditions, we study the bounce and 
singularity formation conditions, junction conditions at 
the boundary of the cloud so as to match it to a suitable 
exterior, and we also consider how the apparent hori- 
zons are affected by the presence of A. One is also inter- 
ested in considering the possibility whether the presence 
of A could restore the cosmic censorship conjecture (see 
e.g. 01 an d references therein), and to investigate how 
the final end state of gravitational collapse is affected by 
A. The present model generalizes the Lemaitre-Tolman- 
Bondi (LTB) dust model studies with A, by intro- 
ducing non-zero pressures in the collapsing cloud. 

The relevant form of Einstein equations, conditions to 
ensure that the collapse develops from a regular initial 
data, and the necessary energy conditions are introduced 
in Section II, together with the details of the tangential 
pressure model. In Section III, we derive the evolution 
of the collapsing matter shells, and explicitly give the 
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conditions when a singularity is formed and when the 
bounce of a particular shell occurs during the collapse 
evolution. The reduction to A ^ dust collapse case, 
when the tangential pressures are put to zero, is also 
demonstrated. For the collapsing solution to be physi- 
cally plausible, it must satisfy certain junction conditions 
at the boundary hypersurface where the interior collaps- 
ing cloud joins with a suitable exterior spacetime. In 
Section IV we study the matching of the collapsing in- 
terior to the exterior Schwarzschild de Sitter or anti-de 
Sitter spacetime, and in Section V we discuss in brief the 
effect of A on the apparent horizons of the fluid model. 
In section VI the nature of the singularity in the tangen- 
tial pressure fluid and dust models is considered when 
A 7^ 0, in terms of its being hidden within a black hole, 
or whether it would be visible to outside observers. The 
final Section VII summarizes some conclusions. 



II. EINSTEIN EQUATIONS, REGULARITY 
AND ENERGY CONDITIONS 

The general spherically symmetric metric in the co- 
moving coordinates (t,r,6,(f>) is given as, 

ds 2 = _ e Mt,r) dt 2 + e 2Ht-r) dr 2 + R 2^ ^2 ^ 

where tf£l 2 = d0 2 + sin 2 6d(j) 2 is the line element on a 
two-sphere. The metric variables e 2v ^ t,r ^ and e 2 '^ t ' r * ) are 
functions of t and r. The e nerg y-momentum tensor for a 
general Type I matter field [12j can be written as, 

Tt = - f r,Tt= Pr ;Tg=T*=po (2) 

The quantities p, p r and pg are the matter density, radial 
pressure and the tangential pressure respectively. It is 
assumed that the matter field satisfies the weak energy 
condition |12|. It means that the energy density as mea- 
sured by any local observer is positive. Then, for any 
timelike vector we have, 

T^V^V > (3) 

This amounts to, 

p + A > 0; p + p r > 0; p + Pe > (4) 

The evolution of the matter cloud is determined by the 
Einstein equations and for the metric these are given 
by (in the units of 8irG = c = 1), 

v'(p + Pr) = 2(pe-Pr)^-p' r (6) 



-2i?' + + R% = (7) 



G - H = 1 - ^ (8) 

where (") and (') represent partial derivatives with re- 
spect to t and r respectively and, 

G(r,t) = e- 2 ^(R') 2 1 H(r, t) = e~ 2v R 2 (9) 

Here F(r, i) is an arbitrary function, and in spherically 
symmetric spacetimes, it has the interpretation of the 
mass function of the collapsing cloud, in the sense that 
this represents the tot al g ravitational mass within a shell 
of comoving radius r [jjj . The boundary of the collaps- 
ing cloud is labelled by the comoving coordinate r n . In 
order to preserve regularity at initial epoch we require 
F(U,0) = 0. 

It can be seen from equation J5J that the density of 
the matter blows up when R = or R' = 0. The 
case R' = corresponds to shell-crossing singularities. 
The shell-cross singularities are generally considered to 
be weak and possibly removable singularities. Hence we 
shall consider here only the shell-focusing singularities 
(taking R' > 0), where the physical radius of all the 
matter shells go to a zero value (R = 0). Let us use the 
scaling independence of the coordinate r to write, 

R(t,r) =rv(t,r) (10) 
where V is the Scale factor. We have, 

v(ti,r) = l; «(t,(r) J r)=0; v<0 (11) 

where U and t s stand for the initial and the singular 
epochs respectively. This means we scale the radial co- 
ordinate r in such a way that at the initial epoch R = r, 
and at the singularity, R = 0. The condition v < 
signifies that we are dealing with gravitational collapse. 
From the point of view of initial data, at the initial epoch 
t — ti, we have five functions of coordinate r given by 
Vq (r) , V'o (f) > Po {r) , Pr ( r ) & Pe i r ) ■ Note that initial data 
are not all mutually independent. To preserve regular- 
ity and smoothness of initial data we must make some 
assumptions about the initial pressures at the regular 
center r = 0. Let the gradients of pressures vanish at the 
center, that is, p' ro {0) = p'g o {0) = 0. The difference be- 
tween radial and tangential pressures at the center should 
also vanish, i.e. p ro (0) — pg (0) = 0. It is seen that we 
have a total of five field equations with seven unknowns, 
P, Pr, Pe, ip, R, and F, giving us the freedom to choose 
two free functions. Selection of these functions, subject 
to the given initial data and weak energy condition, de- 
termines the matter distribution and metric of the space- 
time and thus leads to a particular dynamical evolution 
of the initial data. 

Spherically symmetric collapse models, where the ra- 
dial pressure is taken to be vanishing, but the tangential 
pressure could be non-zero have been studied in some de- 
tail over the past few years 0- The main motivation in 
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the present consideration is to study bounce and singu- 
larity formation conditions for the case when we have a 
non- vanishing A term present, when pressures are also in- 
troduced and allowed to be non-zero within a collapsing 
cloud. One would also like to understand how A affects 
the junction conditions at the boundary where the cloud 
is matched to an exterior spacetime, the horizon forma- 
tion, and the nature of central singularity, when non- zero 
pressures are present within the cloud. 

Taking p r = in Eq. JSJ gives F(t, r) and p in the 
following forms, 



F(t,r) = r 3 M{r) + ^R 3 



3M + r[M, r ) 



v 2 (v + rv') 



(12) 



(13) 



Here A4(r) is an arbitrary function of r subject to the 
energy conditions. There remains the freedom to choose 
one function, since there are six equations with seven 
unknowns. In order to work within the framework of 
a specific class of models, we take v(t, r) in the specific 
form, 



v(t,r) = c(t)+vo(R) 



(14) 



The conditions imposed here, namely that of vanish- 
ing radial pressures, and Eq. i|14fl specifying a form of v 
may be considered to be strong assumptions. However, 
this enables us to make our study in sufficient generality 
with sufficiently rich structure as we shall see below, with 
non-zero pressures introduced into the collapse model. A 
mechanism by which we can have non- vanishing tangen- 
tial pressures is illustrated by the Einstein cluster jlij . 
This is a spherically symmetric cluster of rotating par- 
ticles where the motion of the particles is sustained by 
an angular momentum. This has an average effect and 
creates non-zero tangential stresses within the collapsing 
cloud. 

Also, one can rederive the dust collapse models with a 
non- vanishing cosmological term, by putting vq(FC) = 
in Eq. (|14f> and redefining the comoving time coordinate. 
It is thus clear that the class of models considered gener- 
alize the dust collapse models with a non-vanishing cos- 
mological constant. In general, as v — > 0, p — > oo. Thus 
the density blows up at the singularity R — which will 
be a curvature singularity as expected. Using Eq. (|14H in 
Eq. 0, we have, 



G{t,r) = b{r)e 2 " a{R) 



(15) 



Here b(r) is an arbitrary, C 2 function of r. In correspon- 
dence with the dust models, we can write, 



b(r) = l + r 2 b {r) 



(16) 



where r 2 &o(r) is the energy distribution function for the 
collapsing shells. Using the given initial data and Eq. © 



one can obtain the function vq (R) . Substituting Eq. i|14|) 
in Eq. ||SJ) we get the intrinsic equation of state, 



R 



(17) 



Finally, using equations l|T2* J) .l(T4" j) and (JTBJ in Eq. JSJ), 

we have, 



RR 



a{t)e Uo( ~ R \ (\ + r 2 b )Re 2v ° 



R + r 3 M + ^R 3 



(18) 

Here a(t) is a function of time and by a suitable scaling 
of the time coordinate, we can always make a(t) = 1. Wc 
deal here with the collapse models, and so the negative 
sign is due to the fact that 7? < 0, which represents the 
collapse condition. 



III. SINGULARITY FORMATION AND 
REBOUNCE 

One of our main purposes here is to examine how the 
introduction of a non-vanishing cosmological term mod- 
ifies the collapse dynamics. For example, in the case of 
dust collapse, once the collapse initiates from an initial 
epoch, there cannot be any reversal or a bounce, and the 
gravity forces the cloud to collapse necessarily to a sin- 
gularity. However, this need not be the case when the 
cosmological term is non- vanishing, and we have to re- 
examine the collapse dynamics in order to find how the 
collapse evolves. This we do here for the particular class 
of tangential collapse models as specified above, which 
also generalizes the dust collapse case. 

The evolution of a particular shell may be deduced 
from Eq. ((T51) . Rewriting Eq. (|TB|) in terms of the scale 
factor we have, 



,2v [ v j(r,v) +M + jv 3] 



V(r,v) (19) 



where 



j(r, v) 



(20) 



The right hand side of the equation (|19fl may be thought 
of as an effective potential (V(r,v)) for a shell. The al- 
lowed regions of motion correspond to V{r, v) > 0, as 
v 2 is non-negative, and the dynamics of the shell may be 
studied by finding the turning points. If we start from an 
initially collapsing state [v < 0), we will have a rebounce 
if we get v = 0, before the shell has become singular. 
This can happen when V(r, v) = 0. Hence, to study the 
various evolutions for a particular shell we must analyze 
the roots of the equation V(r, v) = keeping the value 
r to be fixed. It will be seen that the cosmological con- 
stant appearing in the effective potential does play an 
important role in the evolution of a shell. 
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To clarify these ideas, let us consider a smooth initial 
data, where the initial density, pressure, and energy dis- 
tributions are expressed as only even powers of r. Such 
a consideration, that the initial data be smooth, is often 
justified on physical grounds. So we take, 

p(U,r) = Pm+ P2r 2 + Pir A + ■■■ (21) 



pe(ti,r) = p d2 r 2 + P9i r 4 + ■ ■ ■ (22) 



2. If 

b (r) >0; (A + 6p 02 b(r) ) < (27) 

we infer from the sign rule that there is exactly one pos- 
itive root (a(r)). The other two roots are negative or 
complex conjugates. The allowed space of dynamics is 
[0, a]. Thus a > 1 would always ensure a singularity. 
However a < 1 implies an unphysical situation initially, 
where v 2 < 0. 

3. If 



bo(r) = &oo + b 2r 2 + 



(23) 



With the above form of smooth initial data to evolve 
in time using the Einstein equations, we can explicitly 
integrate Eq. © at the initial epoch to get, 



v (R)=pe 2 R 2 



(pe 



P2P0 2 



■R 4 



(24) 



We have neglected here higher order terms in the ex- 
pansion, since at present we want to concentrate on the 
evolution of shells near r = 0. The conditions when 
the treatment is applicable to the whole cloud will be 
discussed later. Eq. (|19|l near the center of the cloud 
(r << r x ) may be written as, 



v 2 = 



(1 + 2pg 2 r 2 v 2 ) [ ( A + 6 P g 2 b(r) ) « 3 + 3fro(r) « + 3A4] 

(25) 

The first factor in V(r, v) is initially positive, because it 
is the | goo | term. As the collapse evolves, the scale factor 
(v) reduces from 1 at the initial epoch to at the time of 
singularity. Hence it clear that the first factor can never 
become zero, and hence does not contribute to a bounce 
of the shell. The main features of the evolution of the 
cloud basically derive from the second factor in V(r,v). 

The second factor in the effective potential expression 
is a cubic equation which in general has three roots. Only 
positive real roots correspond to physical cases. Since the 
coefficient in the second power is zero, we may conclude 
that if all three roots are real then at least one of them has 
to be positive and at least one negative. We observe that 
V(r, 0) = M. > 0. Hence, any region between R = rv — 
and the first positive zero of V(r, v) always becomes 
singular during collapse. The region between the unique 
positive roots is forbidden since in those regions v 2 < 0. 
For a particular shell to bounce it must therefore lie, 
during initial epoch (v = 1), in a region to the right 
of the second positive root. We will now analyze the 
various cases for A ^ in detail and derive the necessary 
conditions. 

1. If 



b (r) > ; ( A + 6pg 2 b(r) ) > 



(26) 



then from the Descartes's rule of signs (See for example 
we see that there are no positive roots. Thus a 
singularity always forms from initial collapse. 



b (r) < ; ( A + 6 P02 b(r) ) < 



(28) 



there is exactly one positive root (f3(r)). Again all shells 
in the allowed dynamical space [0, f3] become singular 
starting from initial collapse. 
4. If 



6o(r) < ; (A + 6 P g 2 b{r)) > 



there are three possibilities, 
4.1 If 

M 2 > 



6g(r 



9(A + 6 P02 b(r)) 



(29) 



(30) 



then there are no positive roots and a singularity is always 
the final outcome of collapse for shells under considera- 
tion. 
4.2 If 



M 2 < -4 



6§(r) 



9(A + 6p 62 b(r)) 



(31) 



then there are two positive roots (71(7") and 72 (r) re- 
spectively). The space of allowed dynamics is [0,71] and 
[72,00). The region (71,72) is forbidden. Shells in the 
[0, 71] region initially, always become singular. But shells 
initially belonging to the region [72,00) will undergo a 
bounce and subsequent expansion starting from initial col- 
lapse. This bounce occurs when their geometric radius 
approaches Rbounce = ^72 • Using the definitions 



Q = 



I -46p(r) 
A + 6pe 2 b(r) 



Cos- 1 



' -9M(r) 2 (A + 6pe 2 b(r)) 
4^ 



(32) 



(33) 



the condition for a particular shell to become singular or 
undergo a bounce may be explicitly written in terms of 
the initial data and A as, 



1 < Q Cosfi ; Singularity 



An 

1 > q Cos{d + — ) ; Bounce 



(34) 



(35) 
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1: The effective potential profile (V(r,v)) for the r = 
with (A + 6pe 2 ) = 7.1978 x 10" 3 , 6 o = -3 x 10" 3 , and 
1.033 x 10 -3 in appropriate units. 



Here g Cos-d is 71 (r) and g Cos{-d + 4^) is 72(7"), which 
are the two roots of the potential function. Note that con- 
trary to the dust models with A, there may be a bounce 
for both positive and negative values of the cosmological 
constant. 
4.3 If 



M" 



b 3 (r) 



%K + Qpg 2 b{r)) 



(36) 



the positive roots are equal. There is no forbidden region, 
and there will be a bounce if 1 > |. Fig. 1. illustrates 
a particular choice of initial data and A which causes a 
bounce in the central (r = 0) shell. 

To analyze the evolution of shells far from the cen- 
ter, in general one has to resort to numerical methods, 
and it is difficult to analytically give a simple expres- 
sion for the singularity or bounce conditions. Neverthe- 
less, the analysis as given here becomes valid for the en- 
tire cloud all the way till boundary when the geomet- 
rical radius of the cloud boundary at the initial epoch 
(R (ti, tv) = rv) is itself small relative to the initial data 
coefficients (ie. p n r™,pg n r™ « 1). Also, if we choose 
the initial data such that the higher coefficients in the 
power series expansion are zero (i.e. Pe n>4 , Pk>4 = 0) 
and p2r 2 ,pg 2 r 2 << 1, then the analysis is again appli- 
cable to the whole cloud. Thus the results derived can 
be considered quite general in these circumstances and 
applicable to the cloud as a whole. 

In this context, if all shells in the collapsing cloud sat- 
isfy Eq. (|3T|) and l|35|> . the complete cloud under- 
goes bounce starting from initial collapse. To avoid shell 
crossings the sufficient condition would be, 



V r e [0,r T ) , 72 (r + S) > 72 (r) 



(37) 



where S is an infinitesimal increment in the comoving ra- 
dius. It is seen that in all the cases discussed, it is not 
A alone, but A + 6pg. 2 b(r) along with b (r) that deter- 
mines the evolution of the shell. This is in contrast to 



the dust models with a non-zero A, where solely the cos- 
mological constant decided the evolution of a shell for a 
given energy function. It is also interesting to note that 
unlike the dust collapse models with A, there could be a 
bounce in the fluid model with vanishing radial pressures 
for both positive and negative values of the cosmological 
constant. This is due to the contribution from the tan- 
gential pressure. 

It can be seen now that one can rederive the known 
bounce conditions in A ^ dust collapse case (Desh- 
ingkar et al 0) as a special case of the consideration 
here, when pg = 0. In that case, 



Av 3 + 3b {r)v + 3M 
3v 



V(r,v) 



(38) 



The study now becomes valid for all shells, without any 
approximation. For example, following similar steps as 
for the fluid model, one obtains for A > 0, 

1. If 

bo (r) > (39) 
the singularity always forms from initial collapse. 

2. If 

60 (r) = (40) 

then again it is found that all shells become singular 
from collapse. 



3. If 

bo (r) < 

there are two scenarios possible. For 

b 3 

M 2 > -4-% 
9A 

all shells become singular. For 



M 2 < -4 



9A 



(41) 



(42) 



(43) 



there are two positive roots (jpi & ^2 respectively) for 
V(r,v). Shells belonging to [0,^1] always become sin- 
gular while those belonging to [^2 > 00 ) undergo bounce 
starting from initial collapse. 



IV. SPACETIME MATCHING 

As we pointed out above, there is a strong physical 
motivation to study and investigate the gravitational col- 
lapse phenomena in a background which is asymptoti- 
cally either a de Sitter or anti-de Sitter metric. For this 
purpose, the collapsing cloud has to be matched at the 
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boundary to a suitable exterior spacetime which has the 
desired properties. 

In the present case, we shall show below that the ex- 
terior vacuum spacetime of the collapsing region may 
be described by the Schwarzschild- de Sitter (SdS) or 
Schwarzschild- Anti de Sitter (SAdS) metric, depending 
on whether the cosmological constant is taken to be pos- 
itive or negative. The collapsing interior cloud which has 
a non-zero tangential pressure is then to be smoothly 
matched to an exterior spacetime in order to generate 
the full spacetime. The necessary and sufficient con- 
ditions to achieve a smooth matching a re g iven by the 
Israel-Darmois j 
shall use below. 



Israel-Darmois junction conditions ( 



lgareg: 

Mm: 



), which we 



Let the interior of the collapsing cloud be described by 
the metric, 

S- : ds 2 _ = -e 2 ^dt 2 + e 2 ^dr 2 + R 2 (t, r)dn 2 (44) 

The exterior vacuum solution can be given as the 
Schwarzschild-de Sitter/ Anti de Sitter spacetime (as de- 
cided by the sign chosen for the cosmological term) as 
given by, 

S + : ds 2 + = -V{TZ)dT 2 + V^iTZ) dTZ 2 + K 2 dn 2 (45) 
where T>(1Z) is given by, 



1 2M \A\TZ 2 



(46) 



The negative sign precedes |A| for a SdS exterior and 
the positive for SAdS exterior. Note that as 1Z — > oo, 
the =F term dominates over the term, and the 
spacetime approaches asymptotically de Sitter metric 
(— ), or the Anti-de Sitter metric (+ 3 ), as the 
case may be. 

Let II denote the boundary hypersurface. The equa- 
tions of the boundary hypersurface considered as an em- 
bedding in the interior or exterior spacetimes are, 

rr:r-iv = o n+ -.n-n^T) = o (47) 

Substituting l|47|) in (|44|l and l|45|l we get the metric on 
the hypersurface as, 

5 n : ds 2 _ = ~e 2u(t ' r ^dt 2 + R 2 (t, r n )dQ 2 (48) 



: ds 2 + = -ViK^d^+V^iK^dKl+Kldri 2 (49) 

The Israel-Darmois conditions to match the interior 
spacetime with exterior require that the first and second 
fundamental forms of the boundary hypersurface match. 
The first fundamental form is given by, 



g^dCdC 



(50) 



where ( parametrizes the hypersurface (£' : t,9,4>). 
The matching of the first fundamental form gives, from 
Eq. |0HJ and Eq. ggjl, 



dr 



3^-^dt 



dr = 



7? 2 



"I 1/2 



dT 



(51) 



(52) 



(53) 



The above three conditions must be satisfied for a smooth 
matching of the collapsing interior to the exterior space- 
time. The next set of conditions will be given by the 
matching of the second fundamental forms. 
The second fundamental form is given by 



Kaud^dC 



where, 



-ntx a 



CSC 



<r 1 (3-y X X^ X 



(54) 



(55) 



is the extrinsic curvature [l8|. The (,) denotes partial dif- 
ferentiation. n a is the normal to the hypersurface which 
is given by, 



f, 



(56) 



where / = is the equation of the boundary hypersur- 
face. Direct calculation gives, 



(0,e 



,0,0) 



(-7T r ,r T ,0,0) 



(57) 



(58) 



The Kgg extrinsic curvatures are calculated as, 

K+= (VKT^ K- e = (RR rr e-^) n (59) 

Now using the fact that the second fundamental forms 
match (i.e. [K^ g — Kqq]^ = 0), we get using equations 
(EJjlEllljlO and (JSnj) after simplification, 

R 2 -2if> n>2 -2v _ -, 2M A| 2 

This is identical to the Cahill and McVittie definitions 
for the mass function (01'E3)j an d hence from Eq. © 
we may by comparison take 



F„ = 2M ± ^-Tll 



(61) 
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This expression suggests that for a smooth match- 
ing of the interior and exterior spacetimes the interior 
mass function at the surface must equal the generalized 
Schwarzschild mass. 

As we can see there is a contribution to the mass func- 
tion from the cosmological constant. A positive cosmo- 
logical constant has an additive contribution and a neg- 
ative cosmological constant has a deductive contribution 
to the mass function F v . The matching of a collapsing 
dust interior with exterior SdS/SAdS would give similar 
results 7]). The K TT and components of the ex- 
trinsic curvature may similarly be calculated. The con- 
dition [Kj T — K~ T ] n = gives no new information since 
the radial pressure is zero. Due to spherical symmetry, 
y^-AA ~ KaAh — gives the same result as Eq. (|6TJ) . 
All the above conditions if5T)) . if5^) . JSBjl and fUlji) must 
be satisfied for smooth matching of spacetimes across the 
boundary. 



V. HORIZONS 

Apparent horizons (H.) are the boundaries of trapped 
regions 01 m the spacetime. We discuss below the effect 
of a non-zero A term on the apparent horizons of the 
tangential pressure fluid collapse models considered here 
briefly. In general, the equation of TL can be written as, 



H: sTR, M R,v=0 
Substituting (g3} in we get, 



2 -74 



R\e 



(62) 



(63) 



From the definition of the mass function Eq. (JHJ), and 
(|63|l we therefore have, 



F 

1 -R=° 



Finally, from O and 

Ti : (3-AR 2 )R = 3r 3 M 



(64) 



(65) 



When A = there is necessarily only one apparent 
horizon given by 



R = r 6 M 



(66) 



which is the Schwarzschild horizon in case we are consid- 
ering that geometry. The same equation also defines hori- 
zon within the collapsing cloud, where R(t, r) is one of 
the metric functions. For the case when A > 0, Eq. I|65|) 
is a cubic equation with at least one positive and one 
negative root. For A > the various cases are given as 
below. 
1. For 



3r 3 M < 



v/A 



(67) 



there are two positive roots for (|65|l and hence there are 
two apparent horizons. These horizons are given by 



R c (r) = -rfCos 
VA 



-Cos- 1 (--r 3 MVX) 
3 2 



Rb(r) = —=Cos 
VA 



f + l -Cos-\-\r 3 M^K) 



(68) 



(69) 



These have been at times called the cosmological, and 
the black hole horizons poj ]. 
2. For 



3r 6 M 



2 

7x 



there is only one positive root for (|65[) . given by 
Rbc{r) 



1 

7x 



This corresponds to a single apparent horizon. 
3. For 



3r 3 M > 



/A 



(70) 



(71) 



(72) 



there are no positive roots and hence there are no appar- 
ent horizons. 

The case l|72fl also shows that the mass of the black 
hole is bounded above by F = and attains the largest 

proper area 47r/A. A general result exists in literature 
|20| showing that in spacetimes with A > and matter 
satisfying the strong energy condition, the area of a black 
hole cannot exceed 47r/A. A detailed treatment of appar- 
ent horizons in the Lemaitre-Tolman-Bondi dust collapse 
models with a non-zero A term is given by Cissoko M. et 
al 0. 

From Eq. I|18f) . the time of apparent horizon formation 
in the fluid model is given by 

tah(r)=t s (r)-t 1 {r) (73) 

where ii(r) is defined as, 

t , , [ Vah{r) sfidv 

h ( r > = / / 

J o Je 4 ' y «vb Q + e 2v ° (v 3 h(rv) + M + jV 3 ) 

(74) 

Also rv ah = R b or R c , h{R) = ( e 2 "°W - l)/R* and t s (r) 
is the time of singularity formation. For shells close to 
the center of the cloud (r << r T ) the expression becomes, 



t a h{r) = t s (r)- 



3>vdv 



/o V ( A + 6 P g 2 ) v 3 + 36 (r) v + 3M 

(75) 

It is observed that the cosmological constant modifies 
the time of formation of horizons and also the time lag 
between horizon formation and singularity formation in 
the fluid model. 
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VI. NATURE OF THE CENTRAL 
SINGULARITY 

The final end state of gravitational collapse and the na- 
ture of the resulting singularity continue to be among the 
most outstanding problems in gravitation theory and rel- 
ativistic astrophysics today. As pointed out earlier, the 
hypothesis that such a collapse leading to a singularity, 
under physically realistic conditions must end in the for- 
mation of a black hole, and that the eventual singularity 
must be hidden below the event horizons of gravity is 
the cosmic censorship conjecture. Despite numerous at- 
tempts, this conjecture as such remains a major unsolved 
problem lying at the foundation of black hole physics to- 
day. 

From such a perspective, we need to examine the na- 
ture of the singularity, in terms of its visibility or other- 
wise for outside observers, when it develops within the 
context of the models considered here. This should tell 
us how the presence of the A term modifies these consid- 
erations, because we already know that in the case of a 
tangential pressure present, but a vanishing cosmological 
constant, both black holes and naked singularities do de- 
velop as final collapse end states depending on the nature 
of the initial data (see e.g. §). 

We have already seen above that the cosmological con- 
stant modifies the time of formation of trapped surfaces. 
If the formation of the horizon precedes the formation of 
the central singularity then the singularity will be nec- 
essarily covered, i.e. it is a black hole. If on the other 
hand, the horizon formation occurs after the singular- 
ity formation, there may be future directed non-spacelike 
geodesies that end in the past at the singularity. Then 
the final end state would be a naked singularity. Thus 
we need to find whether there exist future directed null 
geodesies that end at the singularity in the past. 

Towards analysing this issue, let us define a function 
h(R) as, 

H R ) = W = 2ff(iZ) + 0{R 2 ) (76) 

Using Eq. (|76[) in Eq. (|18|) . we get after simplification, 



rr = - J e 4 "°vb + e 2u ° ( (— + h(rv))v 3 + M ) (7, ) 



Integrating the above equation, we get, 



t(v,r) 



^/vdv 



v yje^vbo + e 2u ° (v 3 (h(rv) + f) + M) 



The time of formation of a shell focussing singularity, for 
a specific shell, is obtained by taking the limits of integra- 
tion in above as (0,1). The shells collapse consecutively, 
one after the other to the center as there are no shell- 
crossings (R! > 0). We are interested in the central shell 
(i.e. the singularity forming at r = 0), since we will see 



that all r > shells are necessarily covered on becoming 
singular. Taylor expanding the above function around 

r = 0, we get, 



t(v,r) = t(v,0) + r 
Let us denote, 



dt(v, r) 



dr 



r=0 



r 2 d 2 t(v,r) 
2! 



d 2 i 



r=0 



X n {v) = 



d n t{v,r) 



dr n 



(79) 



(80) 



The initial data is taken to be smooth (ie. with only 
even powers of r allowed). Due to this choice, the first 
derivative of the functions appearing in the above equa- 
tion vanish at r — 0. Hence we have, 



X x {v) = 



(81) 



Defining B f (r, v) = e iu °vb + e 2v ° ( (f + h)v 3 + M ) we 
may write the Taylor expansion about r — as, 



t(v,r) = t(v,0) 



,- f 1 B f (0,v)y/vdv 



and we have, 



X 2 (v) = 



B f {0,v) 3 / 2 



1 B" f (Q,v)^Jvdv 
B f (0,v) 3 / 2 



(82) 



(83) 



In order to consider the possibility of existence of null 
geodesic families which end at the singularity in the past, 
and to examine the nature of the singularity occurring at 
R = 0, r = in this model, let us consider the outgoing 
null geodesic equation which is given by, 



dt 
dr 



(84) 



We use a method which is similar to that given in |21j . 
The singularity curve is given by v(t s (r),r) = 0, which 
corresponds to R(t s (r),r) = 0. Therefore, if we have any 
future directed outgoing null geodesies terminating in the 
past at the singularity, we must have R — > as t — > t s 
along the same. Now writing Eq. I|84|) explicitly in terms 
of variables (u — r a , R), we have, 



dR 

du 



1 



-(at-l)ft> 



1 - 



r 3 M r AR 2 
R 3 



1 



be 2 " 



Eq. (|85|) is required to be finite and positive, for the exis- 
tence of a naked singularity [21 . In order to get the tan- 
gent to the null geodesic in the (R, u) plane, we choose a 
particular value of a such that the geodesic equation is 
expressed only in terms of (— ). A specific value of a is to 
be chosen which enables us to calculate the proper lim- 
its at the central singularity. In the tangential pressure 
collapse model discussed in the previous section we have 
Ai(0) = 0, and hence we choose a = | so that when the 
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limit r — ► 0, t — > t s is taken we get the value of tangent 
to null geodesic in the (R, u) plane as, 



(86) 



Now note that for any point with r > on the singu- 
larity curve t s (r) we have R — ► whereas, F (interpreted 
as the gravitational mass within the comoving radius r) 
tends to a finite positive value once the energy conditions 
are satisfied. Under the situation, the term F/R diverges 
in the above equation, and all such points on the singu- 
larity curve will be covered as there will be no outgoing 
null geodesies from such points. 

We hence need to examine the central singularity at 
r = 0, R = to determine if it is visible or not. That 
is, we need to determine if there are any solutions exist- 
ing to the outgoing null geodesic equation, which termi- 
nate in the past at the singularity and in future go to a 
distant observer in the spacetime, and if so under what 
conditions these exist. Let xq be the tangent to the null 
geodesies in (R, u) plane, at the central singularity, then 
it is given by, 



,. ,. R dR 
Xo = hm hm — = — — 

t—*t a r^o u du 



t— >i s ;r— >0 



Using Eq. (JBEJl, we get, 



(87) 



(88) 



In the (R, u) plane, the null geodesic equation will be, 

R = x u (89) 

while in the (t, r) plane, the null geodesic equation near 
the singularity will be, 



t-t s (0) 



x r3 



(90) 



It follows that if A^O) > 0, then that implies that 
xo > 0, and we then have radially outgoing null geodesies 
coming out from the singularity, making the central sin- 
gularity locally visible. On the other hand, if #2(0) < 0, 
we will have a black hole solution. We have, however, 
already seen in Eq. Q83[l . that the value of ^(0) entirely 
depends upon the initial data and the cosmological term 
A. Given any A, the initial data can always be chosen 
such that the end state of the collapse would be either a 
naked singularity or a black hole. Hence, it follows that 
for both positive and negative cosmological constants, a 
naked singularity can occur as the final end state of grav- 
itational collapse. 

We noted here earlier that dust collapse models have 
been analyzed in the presence of a cosmological constant. 
The nature of the final singularity in that case has been 
analyzed using the so called 'roots' method (Deshingkar 



S. S. et al 0)- We show below that the different treat- 
ment we have used above to deal with the collapsing 
clouds with pressure included, arrives at similar conclu- 
sions when specialized to the dust case. 

The Einstein equations for dust may be obtained by 
putting pg = in the tangential pressure fluid model 
above. They take the form, 



P(r,t) 



M d (r) 
R 2 R' 



(91) 



(92) 



where f(r) = r 2 bo(r) is the dust energy free function. 
Since we are interested in collapse we must have R <0. 
Then we have from the above equations, 



VRR = -\ M d {r) + f{r)R- 



AR 3 



(93) 



Using the scaling freedom we may write, 

M d {r)= f p(0,r)r 2 dr = p av9 (r)r 3 (94) 
Jo 

Rewriting equation Ij93(l in terms of the scale factor v, 



f (r\v Av 3 

fuv = -\l p av 9 (r) + J -^- + — (95) 

Define ^^} v as L(r, v) which is at least a C 2 function of 
its arguments. Then t s (r,v) is readily calculated from 
to be, 



t s (r,v) 



^/vdv 



v yj ' p av 9(r) +L(v,r) + 



Av 3 
3 



(96) 



The shell focusing singularity R = occurs first for the 
comoving coordinate r = 0. The time of its formation is, 



t.(0) = 



\fvdv 



^/p™f(0) + i(w,0) + 4p 



(97) 



As we did for the fluid model, we Taylor expand t s (r) 
near r — to get, 



t.(r) = i.(0) + rXi(0) + r 2 X 2 (0) + 0{r 3 ) 
where Xi(v) — (assuming smooth initial data), and 



v B d (0,vf/ 2 



with B d (r, v) = p av 9(r) +L(v,r) + ^-. Again, it needs 
to be analyzed whether there exists future directed null 
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geodesies that end at the singularity in the past. Con- 
sider the marginally bound case (/(r) = 0). Then the 
equation of the null geodesic is, 



dr 



In terms of (m, r) it is written as, 



dR 
du 



ar 



a— 1 v 



M d (r) , Ai? 2 



R 



(100) 



(101) 



Choosing a = ^, 

dR 3 ^-Jvv -Rwj 
7 Ar u 



^ + ^) (102) 



Taking the limits as before, the final expression be- 
comes, 



4 3 



1' 



xl = Jpo'MO) 



(103) 



It follows that for both positive and negative values of A, 
there exists initial data that may give X 2 {0) > 0. Thus 
the dust central singularity may be locally visible, even 
when there is a non-zero cosmological constant present. 
It may thus be claimed that the dust central naked sin- 
gularity is not precluded by A. 



VII. CONCLUSIONS 

The final outcome of gravitational collapse is one of 
the most important open problems in gravitation theory, 
and the study of the fluid with vanishing radial pressure, 
and also the dust collapse models here indicate that the 
presence of a non- vanishing A cannot, in any conclusive 
manner, act as a cosmic censor. Although the value of 
the null geodesic tangents are modified by the presence 
of A, there are still cases where the singularity is locally 
visible depending on the initial data. The global visibil- 
ity of such a singularity (null rays from the singularity 
reaching an asymptotic observer) which is locally naked, 
will depend on the overall behaviour of the various func- 
tions concerned which appear in the analysis. But as such 
general relativity has no natural length scale attached to 



it, and so a locally visible singularity may violate the 
cosmic censorship as effectively as a globally visible sin- 
gularity. Moreover, it is seen from various examples that 
if a singularity is locally naked, then one can always make 
it globally visible by suitable choice of the allowed func- 
tions, or by restricting properly the boundary of the cloud 
f[lfH'). We also note that studies pertaining to radiation 
collapse (of a Type II matter field [12j) in spacetimes with 
A, exist in literature Such studies also support the 
conclusion that both collapse endstates, namely a black 
hole or a naked singularity are possible in the presence 
of a non-zero A. 

The currently observed accelerated expansion of the 
universe, and also the AdS/CFT conjecture in string the- 
ories have generated considerable interest in spacetime 
scenarios with a cosmological constant. We studied here 
dynamical collapse of a tangential pressure fluid model 
with non-vanishing A. Various aspects of gravitational 
collapse with A have been discussed in the context of 
the fluid model which include dynamics, junction condi- 
tions, apparent horizons, and nature of the singularity. 
The question of cosmic censorship in the collapse model 
is discussed, and it is also shown that when the pressures 
are put to zero in the model it reduces to the dust case. 

The presence of a finite cosmological constant, it is 
seen, affects the dynamical collapse of the fluid in many 
ways. To summarize, 

1. A plays a role in deciding the dynamical evolution 
of a shell. The conditions for singularity formation and 
bounce are derived. 

2. The matching conditions for the interior and exte- 
rior spacetimes were discussed, assuming that the exte- 
rior is asymptotically AdS or dS. 

3. The apparent horizon conditions and the time of 
horizon formation were studied in brief for A =/= 0. 

4. The effect of A on the nature of the central sin- 
gularity was analyzed. The condition for local visibility 
was derived and it was seen that even for non-vanishing 
A both black hole and naked singularity are possible end 
states of collapse. 
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